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Molecular dynamics simulations and a stochastic Fokker-Planck equation based approach are used to illu¬ 
minate how position-dependent solvent mobility near one or more tracer particle(s) is affected when tracer- 
solvent interactions are rationally modified to affect corresponding solvation structure. For tracers in a dense 
hard-sphere fluid, we compare two types of tracer-solvent interactions: (1) a hard-sphere-like interaction; and 
(2) a soft repulsion extending beyond the hard core designed via statistical mechanical theory to enhance tracer 
mobility at infinite dilution by suppressing coordination-shell structure (Carmer et ah, Soft Matter 2011 , 8, 
4083). For the latter case, we show that the mobility of surrounding solvent particles is also increased by 
addition of the soft repulsive interaction, which helps to rationalize the mechanism underlying the tracer’s 
enhanced diffusivity. However, if multiple tracer surfaces are in closer proximity (as at higher tracer concen¬ 
trations), similar interactions that disrupt local solvation structure instead suppress the position-dependent 
solvent dynamics. 
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I. INTRODUCTION 

Engineering the transport properties of colloidal and 
nanoparticle additives in dispersions is of great funda¬ 
mental interest and has practical implications in a wide 
range of technologies such as drug-delivery mechanisms, 
polymer nanocomposites (PNCs), material fabrication 
techniques, and separations. 1-3 A useful feature of such 
dispersions is that their transport behavior can often be 
systematically varied via tuning the effective inter-species 
interactions, e.g., external electric fields have been used 
to control the dynamics of conducting nanoparticles 4 , sol¬ 
vent pH has been tuned to manipulate the drug-release 
kinetics of nanocapsules 5 , etc. 

One approach to rationally tune transport properties 
is to utilize static-dynamic correlations, such as the em¬ 
pirical and quasi-universal positive correlation between 
excess entropy and long-time particle mobility (e.g., dif¬ 
fusivity) for bulk fluids 6-18 and various inhomogeneous 
fluids 19-28 . Another is through approximate theoreti¬ 
cal approaches, like mode-coupling theory (MCT), which 
connect static structure and dynamic relaxation proper¬ 
ties 29-33 . Within such a framework, one can use statisti¬ 
cal mechanical theories like fundamental measure theory 
(FMT) 34 to connect interparticle interactions and struc¬ 
tural correlations, the latter of which are then approxi¬ 
mately related to the transport properties of the system. 
Using these concepts, one can engineer interparticle in¬ 
teractions for targeted dynamic properties. 

In this vein, previous work using stochastic optimiza¬ 
tion and statistical mechanical theory found that the mo- 
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bility of a single tracer particle in a dense fluid could be 
enhanced (relative to a hard-sphere-like tracer) by adopt¬ 
ing a “flattening” tracer-solvent interaction, which inter¬ 
estingly also increased the effective tracer diameter 35 . In 
particular, the flattening interactions were designed to 
maximize the entropy associated with the tracer’s static 
interparticle correlations, and hence disrupt the coordi¬ 
nation shells of the surrounding solvent particles 35 . How¬ 
ever, several interesting questions emerge from these find¬ 
ings that remain unresolved, including: why should a 
“larger” tracer exhibit faster dynamics? how are the 
dynamics of the solvent particles influenced by tuning 
coordination structure? and what happens to the sol¬ 
vent (and tracer) dynamics if multiple tracers are in close 
proximity, i.e., the average tracer concentration is in¬ 
creased? 

In this paper, we take steps to address these ques¬ 
tions by using a recently introduced steady-state “color” 
reaction-counterdiffusion particle labeling approach 36 to 
calculate position-dependent diffusivities of the solvent 
particles surrounding the tracer additive. This method 
allows us to isolate the affect of the modified tracer on the 
dynamic behavior of the surrounding fluid relative to the 
hard-sphere like tracer, and provides mechanistic insights 
into the underlying phenomenon of how disordering of the 
surrounding fluid affects the tracer-fluid system. 

Furthermore, we explore the impact of tracer concen¬ 
tration on the local coordination-shell dynamics for the 
modified and hard-sphere-like tracer particle, by approx¬ 
imating the solvent-tracer system in a thin-film mor¬ 
phology, where the solvent particles are confined by two 
planar surfaces (representing the tracer particle surface, 
Fig. 1). This model is motivated by experimental stud¬ 
ies that have also shown the utility of approximating the 
bulk material properties of high-additive composites via 
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FIG. 1. (color online). Schematics illustrating the two system 
geometries under consideration, (a) Tracer particle (t) at in¬ 
finite dilution surrounded by a 3D volume of dense solvent (s) 
particles, where we consider solvent motions along the radial 
r-coordinate relative to the tracer center of mass, (b) Thin 
film approximation of systems with higher tracer concentra¬ 
tion p t , where solvent particles are confined between walls 
(w) representing tracer surfaces. The walls are separated by 
thickness H in the 2 -direction, and periodic boundaries are 
employed in the x- and ^-directions. Particles in (a) and (b) 
are drawn on the same scale. 


examination of thin films having similar surface proper¬ 
ties. 37,38 . 


II. COMPUTATIONAL METHODS 


pair potential 40 : <p ss (r) = <p WCA {r,c r s ,e) = 4e([<r s /r] 48 — 
[cr s /r] 24 ) + e for r < 2 1//24 <r s and <p BB (r) = 0 for r > 
2 1 / 24 cr s , where r is the interparticle separation. The 
tracer-solvent pair interaction is defined by: <pts(r) = 
ip WCA ( rjC r ts)e ) for r < r 0 and ip ts (r) = ip WCA (r, a tB , e) + 
<Pts,o(»") for r > r 0 , where cr ts = (cr t +cr s )/2. Here, tpts,o( r ) 
is a contribution outside the hard core that can be tuned 
to affect the solvent density profiles around the tracer, 
and ro is chosen such that the tuning procedure does not 
alter the range of the hard-core exclusion interaction. 35 

For the thin films in Fig. 1(b), we simulate N s = 
4000 solvent particles confined between between two flat 
walls (w) in the 2 -direction separated by a thickness 
H , and are periodic in the x- and y-directions. The 
wall-solvent (or tracer-solvent) interaction is defined by 
V?ws(z) = + WCA (* + 0.5, er ss ,e) for z < z 0 and <py, B (z) = 
ip WCA (z + 0.5, cr sa ,e) -(- Ws,o(z) for z > z 0 - This term 
is equivalent to the solvent-solvent hard-core interaction 
combined with an additional variable </Vs,o(z) that can 
be tuned to affect the solvent density variations across 
the film. Here, we define <^ W s(z) as a function of (2 + 0.5) 
such that solvent centers can access positions ranging 
from approximately 0.5 < z < ( H — 0.5), as is the case 
for true hard spheres situated between flat hard walls. 

We are primarily interested in measuring how solvent 
motions near the tracer and wall surfaces depend upon 
the type of density profiles allowed in the near-surface sol¬ 
vation layers. Notably, we are probing diffusive solvent 
displacements along paths that are inhomogeneous-i.e., 
motions along the r- and 2 -coordinates that are subject 
to non-isotropic solvent density fields-such that we can¬ 
not estimate the corresponding position-dependent sol¬ 
vent diffusivities from particle displacements using the 
typical Einstein relation 41 . (However, the Einstein rela¬ 
tion is applied to obtain, e.g., isotropic tracer diffusivity 
D t = (Ar 2 )/6A t, where (Ar 2 ) is the tracer mean squared 
displacement in the x, y, and 2 directions over lag times 
At exceeding the timescales of ballistic motions.) 

Position-dependent particle diffusivities in dense inho¬ 
mogeneous fluids are accurately described by the Fokker- 
Planck (FP) equation 36,42-44 . For the single-tracer case 
in Fig. 1(a), the FP equation describing solvent displace¬ 
ments along the r-coordinate of the (reference) tracer 
particle is 42 


For both the single-tracer and thin film systems shown 
in Fig. 1, we generate particle trajectories from 3D molec¬ 
ular dynamics (MD) simulations in the canonical ensem¬ 
ble, where trajectories are evolved using the velocity- 
Verlet method 39 . We use a time-step of O.OOlcr s \Jm/e, 
where <j B , m, and e are the characteristic solvent (s) diam¬ 
eter, mass, and energy scales, respectively. The tempera¬ 
ture is constrained to T = e/k b (where ks is Boltzmann’s 
constant) using the Nose-Hoover thermostat. 

For the case in Fig. 1(a), we incorporate a single tracer 
(t) with diameter a t into a fluid comprised of N s = 8000 
hard-sphere (HS) solvent particles. Solvent-solvent HS 
interactions are approximated by the following continu¬ 
ous, steeply-repulsive Weeks-Chandler-Andersen (WCA) 


dG 

~dt 


1_d_ 

r 2 dr 




(1) 


which contains position-dependent diffusivities in the re¬ 
direction D* (r). Here, G(r,t 0 + At\r',t 0 ) is the Marko¬ 
vian propagator describing temporal single-particle dis¬ 
placements given a non-uniform potential of mean force 
(PMF). For the solvent surrounding a single tracer, this 
PMF is given by V^ B (r) = - ln{[p t s/Ats,avg](?*)} + C, 
where pts{r) is the solvent density at some distance r; 
[pts /pts,avg] {p ) is the partial radial distribution function 
(i.e., <7ts(0)> an d C i s an arbitrary constant. 

To extract D S I (r) profiles from particle trajectories, 
we use a color reaction-counterdiffusion treatment of the 
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FIG. 2. (color online), (a) Radial density profiles of solvent 
particles [p s /p s , avg ](r) surrounding tracer particles of diame¬ 
ter at/a a = 2.0 for the two different tracer-fluid interactions 
<pts(r) shown in inset (b). Red solid ‘hard-sphere’ profiles cor¬ 
respond to bare WCA potential for <p ts (r), while blue dashed 
‘flattening’ profiles correspond to a ipts(r) optimized to flat¬ 
ten [ps/p s ,avg](?’) and thus maximize the contribution of the 
tracer particle’s correlations to the two-body excess entropy 
associated with tracer-solvent static correlations, (c) Cor¬ 
responding position-dependent diffusivities of fluid particles 
D s>r (r) along the radial direction. 

steady-state form of eq. 1 (i.e., dG/dt = 0) 36 . We 
then rescale the D StT (r) profiles that account for the rel¬ 
ative mobilities of the the tracer particles themselves, 
where D SjI (r) = D* T (r) [D“ k /(Dg lk + J9 t )], A,r( r ) are 
the non-normalized local diffusivities obtained from eq. 1, 
_D klk is the solvent bulk diffusivity in the absence of any 
tracer, and D t is the long-time tracer diffusivity. Thus, 
As,rM approaches _D klk as r —>■ oo independent of the 
type of tracer-solvent interaction. Implementation de¬ 
tails and the analogous expressions for local solvent dif¬ 
fusivities along the z-direction _D sz (z) for the thin film 
systems can be found in a previous publication 36 . 


III. RESULTS & DISCUSSION 

We begin our discussion by considering Fig. 2, where 
we demonstrate how solvent structure and dynamics 
around infinitely dilute tracer particles are affected by 
two different types of tracer-solvent interactions <p ts 
shown in Fig. 2(b): (1) a hard-sphere-like steep WCA 
repulsion (i.e., the contribution outside the hard core 
Ptfiir) = 0 for r > ro); and (2) the same core repulsion 
combined with a softer long-range repulsion (resembling 
Yukawa screened electrostatic interactions) to flatten 


(i.e., eliminate) solvent coordination shells around the 
tracer. As shown previously 35 for a single tracer, adopt¬ 
ing such flattening potentials minimizes tracer-solvent 
structural pair correlations which increases the cor¬ 
responding two-body excess entropy S 2 ,ts relative to the 
hard-sphere-like tracer case, as evident from its defini¬ 
tion: 

f‘00 

S2,t s /A = -y J {gt s (r)\ng ts (r) - g ts {r) + l}dr (2) 

In turn, the long-time tracer diffusivity D t can be en¬ 
hanced by a factor of up to two or more compared to the 
hard-sphere-like WCA case depending on the size ratio of 
the tracer and solvent particles at/cr s . 35 This is-at first 
glance-counterintuitive because the tracers with the flat¬ 
tening potentials have larger apparent diameters, which 
one might naively expect to depress diffusive mobility. 

In Fig. 2(a) and 2(c), we show the position-dependent 
solvent densities [p s /As,avg]M (be., fftsfj)) and diffusiv¬ 
ities A,rM around these two tracer particles of diam¬ 
eter C7 t /a s = 2.0. The average tracer diffusivities from 
both cases were related as: D^ at /D^ VGA ~ 2 (see Fig. 
3 in ref. 35 ). Here, it is evident that using the flattening 
V’ts.oM thoroughly destroys solvent coordination shells 
while simultaneously enhancing local solvent diffusivities 
D s , r (r) at virtually all distances at and near the tracer 
surface, which helps to rationalize the corresponding en¬ 
hancement in tracer diffusivity. The solvent diffusivity 
measurements for the two cases reflect known excess en¬ 
tropy scalings for transport coefficients and support the 
notion that tracer and solvent dynamics are coupled, i.e., 
the shift in solvent dynamics drives (even non-intuitive) 
trends in relative tracer mobility. 

While the discussion above shows that destroying coor¬ 
dination shells enhances tracer and solvent dynamics for 
infinitely dilute tracers, a natural question is whether the 
use of this strategy would have the same implications for 
systems with finite tracer concentrations p t > 0. To ad¬ 
dress this, we consider how solvent dynamics are affected 
in the thin film systems illustrated in Fig. 1(b), where the 
interwall separation H approximately scales with tracer 
concentration according to H oc p t ' .In this way, we ex¬ 
amine whether destroying density variations categorically 
enhances nearby solvent diffusivity-and tracer diffusivity, 
due to their coupling-regardless of the average distance 
between proximal “tracer” surfaces. 45 

In Fig. 3, we show the position-dependent solvent den¬ 
sity p s (z) and diffusivity A,z(~) profiles associated with 
wall-solvent interactions <p ws (z) that are either steeply 
repulsive (i.e., ip w fi(z) = 0) or have been augmented 
with <p w ,o(~) potentials derived from Fundamental mea¬ 
sure theory (FMT) 22,34 to flatten the solvent density pro¬ 
files p s (z) across different film thicknesses H. Based 
on the top panels of Fig. 3, it is evident that ramp- 
like <p Wj o (z) repulsions spanning the characteristic parti¬ 
cle lengthscale a s are effective at destroying near-surface 
density variations. Such flattening interactions also re¬ 
semble the ipt,o(z) flattening interactions for spherical 
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FIG. 3. (color online) Results for thin film systems: (top) 
wall-solvent interactions /3ip wa (z); (middle) solvent density 
profiles p s (z) normalized by average density p B}avg ; and 
(bottom) position-dependent diffusivities of solvent particles 
D SiZ (z) in the 2 -direction. Left panels correspond to film 
thickness H = 5 and right panels to H = 15, where in both 
cases the average solvent packing fraction is <j> = 0.35 (based 
on total channel thickness H). Red solid ‘hard-sphere’ pro¬ 
files correspond to the WCA potential for fhp wa (z), while blue 
dashed ‘flattening’ profiles correspond to /?<^ ws ( 2 ) potentials 
derived via FMT to flatten p B (z) across H. 
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FIG. 4. (color online) (a) Average fluid diffusivities in the 
2 -direction D a ff(z) for various thin film thicknesses H/o a , 
where D a ff(z) is calculated based on 2 -positions near the 
confining surfaces (see text), (b) Average diffusivities in the 
unconfined directions D xy for fluid particles within the same 
regions, (c) Average two-body excess entropies quanti¬ 
fying static correlations between solvent particles within the 
same regions. Lines connecting symbols are a guide to the 
eye. 


tracers. By examining the middle- and bottom-left pan¬ 
els of Fig. 3, one can see qualitatively different dynamic 
responses to flattening p s (z) for various H: the diffusiv¬ 
ities D s , z (z) are uniformly depressed for the highly con¬ 
fined H = 5 case while they are instead enhanced near 
the walls for H = 15. The latter observation is expected 
as H —> oo qualitatively corresponds to the single-tracer 
limit p t —>• 0. 

We can understand these opposing results in a gen¬ 
eral way by considering particle packing effects within 
thin films. In the case of unmodified hard-wall-like WCA 
boundaries, particles tend to accumulate near the walls 
to minimize excluded volume. Given that the external 
potentials capable of flattening p s {z) are soft repulsions 
that “push” these particles toward the center of the film, 
it is perhaps unsurprising that for very thin films, there is 
insufficient space for these particles to redistribute them¬ 
selves in a way that allows for efficient packing and cor¬ 
responding diffusive motions 22 . In contrast, as the thick¬ 
ness of the thin film increases, particles near the walls 
constitute a smaller fraction of the total fluid popula¬ 
tion. Thus, their redistribution away from the walls is 
accomodated more readily, which should not uncondi¬ 
tionally frustrate near-surface diffusion. 

In Fig. 4, we demonstrate that this dichotomy in dif¬ 
fusive responses to flattening p s (z) at various H is sys¬ 


tematic and can be rationalized by considering the cor¬ 
responding changes in the static solvent-solvent correla¬ 
tions, as quantified by excess entropy. To make this anal¬ 
ysis consistent, we calculate average solvent diffusivities 
both perpendicular D and lateral to the hard- 

sphere-like and flattening potentials for particles located 
in regions near the thin film walls. These regions span 
0 < z < d ws and (H — d ws ) < z < H , where d ws = 2.5a s 
for all values of H (this d ws corresponds to complete cov¬ 
erage for the H/a s = 5 film). We also calculate average 
two-body excess entropies sj^ 8 using the same particle 
trajectories within this region. 

Given D SjZ (z) diffusivity profiles, average diffusivities 
in the near-surface regions are calculated via 

D Z e = \ f D sA z )p( z ) dz / [ P( z ) dz ( 3 ) 

i=o,i Zi Zi 

where zq = 0, z' 0 = d ws , z\ = H — d ws , and z[ = H. 
We calculate lateral diffusivities via the Einstein relation 
-^s,xy = (Ar 2 )/4A t using mean squared displacements 
(Ar 2 ) in the periodic x and y directions. Displacements 
are aggregated from any particles that are located within 
d ws of the nearest film boundary at time t up until the 
maximum time lag Af max = t' — t 1 where t 1 is the time 
at which the particles exit the near-wall region. 46 We 
obtain values via an integral analogous to eq. 3 over 
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position-dependent two-body excess entropy S 2 (z), where 
S 2 (z) profiles are calculated using a recast form of eq. 2 
in the Cartesian geometry 44 . 

In Fig. 4(a-b), we compare average solvent diffusivities 
measured near surfaces with ‘hard-sphere-like’ and ‘flat¬ 
tening’ ip ws (z) interactions for various H , where it is evi¬ 
dent that the effect of flattening density variations upon 
solvent mobility is qualitatively dependent upon inter¬ 
surface proximity. For Ha s < 10, average dynamics in 
all directions are slowed upon flattening the density pro¬ 
file, while for H/cr s > 10, dynamics are instead enhanced 
(consistent with the tracer case). This crossover can be 
unified with the single-tracer case by considering the 
curves shown in Fig. 4(c), which illustrate that the rel¬ 
ative slowdown for H/a s < 10 reflects the more gen¬ 
eral positive correlation between particle mobility and 
multi-body excess entropy. Thus, eliminating one-body 
density variations has a non-trivial FZ-dependent effect 
upon two-body static correlations, where the latter are 
more meaningfully correlated with dynamics. (For an¬ 
other pronounced example of such effects, see Goel et 
al. 22 .) 

Crucially, the trends in Fig. 4 imply that tracer- 
solvent interactions <pts(r) designed at infinite dilution 
to enhance tracer (and solvent) mobility may not gen¬ 
erally have the same qualitative impact at sufficiently 
high tracer concentrations (i.e., sufficiently thin inter¬ 
tracer solvent regions). Using the crossover film thick¬ 
ness H* /a s ~ 10 from Fig. 4, we can obtain an order-of- 
magnitude estimate for the limiting tracer concentration 
pi beyond which the single-tracer physics might be ex¬ 
pected break down due to solvent packing effects: our 
most conservative (i.e., biased towards a greater value of 
p*) calculations indicate that pi « 0(0.01). 47 That pi 
is so low lends a cautionary note in terms of deploying 
tracer-solvent interactions designed in the dilute limit at 
any significantly higher concentrations. 


IV. CONCLUSIONS 

Using recently introduced techniques for characteriz¬ 
ing the position-dependent dynamics of inhomogeneous 
fluids, we have shown that eliminating the coordina¬ 
tion structure of an infinitely dilute tracer (additive) 
in bulk solvent to decrease static correlations (increase 
excess entropy)-which is achieved by rationally tuning 
the tracer-solvent pair interaction-enhances the diffusive 
mobilities of both the tracer and the surrounding sol¬ 
vent particles. However, upon incorporating similarly 
tuned interactions into thin films of solvent particles, 
which approximate systems of higher tracer concentra¬ 
tion, we find that eliminating one-body solvation struc¬ 
ture for film thicknesses (i.e., tracer-tracer distances) 
smaller than several solvent particle diameters decreases 
excess entropy and suppresses solvent mobility. This ob¬ 
servation nicely explains the results of two previous stud¬ 
ies where tuning interactions to increase excess entropy 


had opposite effects on the dynamics of solvent surround¬ 
ing a tracer particle 35 or trapped between interacting sur¬ 
faces 22 . Due to the apparent coupling of tracer and sol¬ 
vent dynamics, this is suggestive that tracer-solvent in¬ 
teractions designed at dilute conditions could have qual¬ 
itatively different impacts upon system dynamics at suf¬ 
ficiently high tracer concentrations, where the critical 
loadings separating these regimes are likely very low. 
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